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I. INTRODUCTION 



When the temperature T of the quark-gluon plasma is much greater than the QCD scale 
parameter Aqcd, the hard modes, i.e. those with momenta of the order of T or larger, are 
weakly interacting and they can be described within perturbative QCD |]]. The dynamics 
of the soft sector, however, remains non-perturbative even at arbitrarily large temperature 
0], as signalled by severe infrared divergences ||. Then, one has to refer to effective theories 
to get an insight into the soft mode dynamics. Such theories, see e.g. can be derived 

from QCD by integrating out the hard modes, but constructing them is by far not a simple 
task. Consequently, one often relies on more or less heuristic approaches, usually exploiting 
a semi-classical or classical field approximation because the occupation numbers of the soft 
gluonic modes are large. 

A very natural effective approach is provided by the kinetic theory, where the hard modes 
are treated as (quasi-)particles while the soft gluonic ones contribute to the chromodynamic 
mean field. The transport theory has been formulated in two versions. The first one 
treats the color degrees of freedom as a classical continuous variable which, as position or 
momentum, evolves in time. A starting point of the theory are the Wong equations |J, 
which describe a classical particle that interacts with the chromodynamic field due to the 
color charge. Then, one immediately gets the Liouville and the transport equations [|7| of 
a many-body quark-gluon system. The physical content of the theory is rather transparent 
and numerous results, for example transport coefficients, can be easily obtained. Even the 
simplest collisionless transport equations, where the dissipation phenomena are neglected, 
provides a surprisingly rich dynamics. The transport theory with the classical color became 
really reliable when the theory was found H to reproduce the QCD hard-thermal-loop 
dynamics P,p|,|T0|. It was further established (Tl|] that the theory supplemented by the 



collision terms, obtained integrating out soft fluctuations around the mean fields, agrees 
with the QCD effective approaches f5,12fl. The relationship between the transport theory 
with classical color and QCD can be found through the study of the quantum path integral 
within a saddle-point approximation fll3| . 

In the second version of the QCD transport theory |]7|,^4|, the color charges are repre- 
sented, in full accordance with QCD, by a matrix structure of the distribution function. 
The Vlasov transport equation of quarks was derived [l5j directly from QCD, by analyzing 
the motion of quantum quarks in the classical chromodynamic field. The gluon transport 
equation was found |16|JT7|| , by splitting the gluon field into the mean field and the contri- 
bution representing the particle excitations. It was further observed [0] that the quark and 
gluon transport equations are formally identical when the first one is written in the funda- 
mental representation and the second one in the adjoint representation. Then, the quark 
and gluon distribution functions are N c x N c and (N% — 1) x (jVj? — 1) matrices, respectively, 
for the SU(N C ) gauge group. The early development of the quark-gluon kinetic theory was 
summarized at ref. |19| . 

In quasi-equilibrium, the matrix transport theory was proved |2(J to be fully equivalent 
to the QCD hard-loop approach 0|],TI]]. The kinetic approach, which can be treated as 
a local representation of the non-local hard-loop action, is particularly useful to study the 
collective excitations of the quark-gluon plasma, see Gil] for a review. More recently, the 
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quasi-equilibrium kinetic equations have been derived beyond the collisionless limit [22] and 
the QCD effective theories || have again been correctly reproduced. 

A natural question that arises is where the agreement between the kinetic theory, either 
with the classical color or in the matrix form, and the finite temperature diagrammatic 
approach breaks down. Surprisingly enough, the effective action of the static fields provided 
by the kinetic equations agrees with that obtained within perturbative QCD even at the g 3 
order , where an operator responsible for C-odd processes appears for systems with finite 
baryon density. However, in a subsequent study p3j the kinetic theory with the classical 
color has been found to reproduce the g 4 contribution to the effective action only in the limit 
of high-dimensional color representations. Thus, the limitations of the classical approach 
have been explicitly determined. 

The aim of this paper is to clarify whether the difficulties faced by the classical color 
transport theory can be overcome when the matrix formulation is used. We explore how 
far the limits of the non-Abelian kinetic approach can be extended. For this purpose, we 
look for the solutions of the transport equations for quarks and gluons interacting with a 
chromodynamic mean field. The system is assumed to be translation invariant in one or 
more space-time directions. Thus, our considerations hold, in particular, for static and for 
homogeneous systems. We first present exact solutions when the vector gauge fields in the 
direction of the translation invariance commute with their covariant derivatives. Then, we 
use perturbation theory to solve the transport equation when the commutation condition is 
relaxed. Once the solutions are known, we derive the corresponding color current, and then 
the associated effective action. In the case of thermodynamic equilibrium, our results agree 
with those obtained by computing the one-loop effective action of QCD in the presence 
of a constant background field, which also corresponds to the effective potential of the 
dimensionally reduced theory [27]. However, for static and non-constant background fields 
such that the commutation condition is not satisfied, we find additional non-local operators 



that correct the effective potential of H27| , in what seems to be a discrepancy between the 
two approaches. 

The paper is organized as follows. In Sec. [H| we review the transport theory approach 
we use in this paper. Exact solutions of the transport equations for translation invariant 
systems are discussed in Sec. |T| while in Sec. [TV] we solve the equations perturbatively. The 
effective actions for static and for homogeneous systems close to equilibrium are derived in 
Sec. [V] and we conclude our considerations in Sec. [V]]. The evaluation of some momentum 
integrals is left for Appendix A and we collect in Appendix B some formulas of the traces of 
SU (N c ) generators in the adjoint representation. Finally, in Appendix C we briefly discuss 
the one-loop effective action P,^5|,|26| which is compared to our results. 



II. TRANSPORT EQUATIONS 

In this section we briefly review the transport theory of quarks and gluons [19|jll| . While 
we will restrict the discussion to QCD, with Nf massless quarks and antiquarks carrying 
color in the fundamental representation and gluons in the adjoint, the results could easily 
be generalized to a different non-Abelian theory with different field content. 
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The distribution function of (anti-)quarks Q(p, x) (Q(p, x)) is a hermitian N c x N c matrix 
in color space (for a SU(N C ) color group); x denotes the space-time quark coordinate and 
p its momentum, which is not constrained by the mass-shell condition. The spin of quarks 
and gluons is taken into account as an internal degree of freedom. The distribution function 
transforms under a local gauge transformation M as 

Q{p,x) -> M(x)Q(p,x)M\x) . (2.1) 

Here and most cases below, the color indices are suppressed. The distribution function of 
hard gluons is a hermitian (iV| — 1) x (AT? — 1) matrix, which transforms as 

G(p, x) -> M(x)G(p, x)M\x) , (2.2) 

where 

M ab (x) = Tr[r a M{x)nM\x)\ , 

with r a , a = l,..., N% — 1 being the SU(N C ) group generators in the fundamental represen- 
tation with Tr(r a T b ) = |5 a6 . 

In a collisionless limit, the distribution functions of quarks and gluons satisfy the trans- 
port equations: 

p^D,Q(p,x) + ^p^F, u (x),^^^ = 0, (2.3a) 

P^D,Q( P ,x)-^p^If, u (x),^^-^ = 0, (2.3b) 

^ M G(p,x) + |^|^(x),^^| = 0, (2.3c) 

where g is the QCD coupling constant, {...,...} denotes the anticommutator; the covariant 
derivatives _D M and act as 

£> M = <9 M - iglA^x), ... ] , £> M = <9 M - ig\A^x) } ... ] , 

and A^ being four-potentials in the fundamental and adjoint representations, respec- 
tively: 

A»{z) = A£{x)r a , M^x) = -iUA^) , 

and fabc are the structure constants of the SU(N C ) group. Since the generators of SU(N C ) in 
the adjoint representation are given by (T a ) bc = —if a bc, one can also write A 1 * = A£T a . The 
stress tensor in the fundamental representation is F^ v = d fl A u — d u A fl — ig[A^,A u ], while 
T^y denotes the field strength tensor in the adjoint representation. 

Sometimes it is convenient to project the matrix equations ( |2.3| ) into their colorless and 
colored components. For the quark distribution function we define 

Q(p, x) = q(p, x) + q a (p, x)r a . (2.4) 

Then, we can deduce from Eq. ( |2.3a|) a set of coupled equations for the different components 
of Q defined by Eq. (|2.4|). More precisely, we find 



p»d»q(p,x) + A^F^ix)^^- = , (2.5a) 
p^Vip, x) + 9 - d^ u (x)?^l + 9P"F^x) d -^L = , (2.5b) 



where d abc are the totally symmetric structure constants of SU(N C ) and D^f = d^5 ac + 
gjabcj^b Similar equations can be written for the antiquark and gluon distribution functions. 
Eqs. Q2.5| ) reflect in a very clear way that transport phenomena of colorless and colored 
fluctuations are coupled beyond the lowest order in the gauge coupling constant. Eqs. (|2.5| ) 
might be very useful when collisions are also taken into account. In this study, however, we 



find more convenient to work with the matrix equations (2.3). 

Once the solution of the transport equations is known, we can obtain the color current 
associated to the plasma constituents. The color current is expressed in the fundamental 
representation as 

f{x) = "f / dPp»[Q(p,x) - Q(p,x) - ^Tv[Q(p,x) - Q(p,x)} + 2ir a f abc G bc (p,x)] , 

(2.6) 

so that j%(x) = 2Tr(r a j M (x)), and the momentum measure 

dP = -0- 3 2Q(p o )6(p 2 ) (2.7) 

takes into account the mass-shell condition p = |p|. Throughout the paper, we neglect the 
quark masses, although those might easily be taken into account by modifying the mass- 
shell constraint in the momentum measure. A sum over helicities, two per particle, and over 
quark flavors Nf is understood in Eq. (|2.6|) , even though it is not explicitly written down. 

In the transport theory framework one can consider two different physical situations: 
1) the gauge fields entering into the transport equations ( |2.3| ) are external, not due to the 
plasma constituents; 2) the gauge fields can be generated self-consistently by the quarks 
and gluons. In the latter case, one also has to solve the Yang-Mills equation 

D,F^(x)=f(x), (2.8) 



where the color current is given by Eq. ( p.6|) . 

The color current can be derived from an effective action added to the Yang-Mills one. 
By means of the relation 

it - • (2-9) 

where S = J d A x £ , one can obtain the effective action, up to an integration constant, from 
the knowledge of the color current. In the remaining part of this article we will use this 
approach to obtain the effective action in different physical situations. 
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III. EXACT SOLUTIONS 



Finding exact solutions of the transport equations Q2.3Q is in general a difficult task. 
However, it is possible to find such solutions under some restrictive conditions. Here we 
consider a system where both the vector gauge field and the distribution functions are 
invariant under the space-time translation(s), i.e. 

d ai A»{x) = 0, /i = 0,l,2,3, (3.1) 

and 

d ai Q(p, x) = d ai Q(p, x) = d a .G(p, x) = , (3.2) 

for a fixed aj, where oti can involve more than one Lorentz index. For example, if a,i = the 
system is static while for a, = 1, 2, 3 the gauge field and the distribution functions depend 
only on time. The condition Eq. (|3.1|) is a choice of gauge that, as we will show below, allows 



one to find solutions to the transport equations that respect the translation invariance of 
the system. 

To solve Eqs. ( |2.3| ) for translation invariant systems along the direction x ai , we will take 
into account two known facts. First, in an electromagnetic plasma, an exact solution of 
the corresponding transport equation is given by any function of the canonical momentum 



p ai — eA a -(x), if d ai A^ = p8[ . Second, the non-Abelian transport equations of particles 
carrying a classical color charge I a for these translation invariant systems are also solved 
by any function of the canonical momentum p a . — gA^.(x)I a , if d^A^ = f2"3 |. Because 



the first case represents the Abelian limit of the transport equations (|2.3i), while the second 
one corresponds to the limit of high- dimensional color representations^, one expects that 
the solutions of Eqs. fl2.3|) are of the form: 

Q(p,x) = fba, ~ gAa.ix)) = £ y -J^A ai (x) A a2 (x) ■ ■ ■ A an (x) a JKPai > , (3.3a) 

Qip, x) = fipa, + gA ai = £ ^jA ai (x) A a2 (x) ■ ■ ■ A an (x) , (3.3b) 



G(p,x) = f 9 ( Pat - gA ai (x)) = E K -^A ai {x)A a ,{x) ■ ■■A otn {x) , ^ /g(Pa< l , (3.3c) 

where it is understood that a sum is taken over the repeated indices. The functions /, / and 
f g are, in principle, arbitrary but they can be fixed by additional considerations. 

Let us note that the distribution functions given by Eqs. (|3.3|) transform covariantly, i.e. 
according to Eqs. Q2-lp and ( |2.2| ), even though the potentials A^, A^, in general, do not. 
Indeed, 



tWhile we are considering here the transport equations only for particles carrying color in the 
fundamental and adjoint representations, the equations are expected to have the same structure 
in any other non-Abelian representation. 
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A a Hx) -> M(x)A a Hx)MHx) - -(d ai M(x))M ] (x) . 

9 



(3.4) 



However, the second term in the r.h.s of Eq. ( |3.4|) , which transforms non-covariantly, is 
eliminated because of the condition (|3.1|) . 

Let us check under which conditions the ansatz Eq. (|3.3| ) solves the transport equations. 
We first note that Eq. ( |3.3aD is totally symmetric under the exchange of indices a±, ■ ■ ■ , a n . 
Inserting Eq. fl3.3a| ) into the transport equation (|2.3aj) , one finds 



oo / \n n— 1 

jfDpQfax) = E 4* ■ ■ ■ A ag (D^A as+1 ) ■ .-A ar , 



n=0 



III 



s=0 



dp ai dp a2 ...dp 



(3.5) 



9 L 

2 y 



dQ(p, x) 



:-g) 



n+l 



2n\ 



\^D^A ai , A ai ■ ■ ■ A an | 



d n+1 f(p 



d Vv J n=0 
where we have used the following property of a commutator 



dp ai dp ai . . . dp an 



(3.6) 



n-l 



[X,Y n ] = Y,Y S [X,Y]Y 



n—s—l 



s=0 



to derive Eq. (|3.5|). We have also taken into account that F^ a . = D^A a . because of Eq. (|3TT|). 

It is easy to observe that the terms in Eq. ( |3.5|) and Eq. (|3.6|) that correspond to the 
same order of the derivative of / cancel each other exactly if 



[D^Aa- , A a . 



0, 



0,1,2,3 



(3.7) 



where A a . is also in the direction of the translation invariance. The same condition is 
obtained for the antiquark distribution function, while for gluons one finds that Eq. ( |3.3cj ) 
is an exact solution of Eq. (|2.3c|) if 



[2VUA*J=0, /i = 0,l,2,3. (3.8) 

However, it is easy to prove that Eq. ( ^8|) is automatically satisfied if Eq. (|3.7| ) holds. 

For static systems [a = 0), Eq. ( |3.7| ) reduces to the commutation relation between the 
color electric field and A . In a more general situation, the condition Eq. ( |3.7|) simplifies 
the non-Abelian field dynamics in the direction of the translation invariance. Note that 
the commutation condition is trivially satisfied in the Abelian limit. If is a generator of 
a representation R of SU(N C ), then, after a normalization of these generators, one would 
get [TjJ, Tjy — > for high-dimensional representations. Consequently, the commutation 
condition would also be satisfied. This explains how to reconcile the matrix results with 
those obtained with the non-Abelian transport equations for classical color. 

Once the solution for the quark, antiquark and gluon distribution functions are known, 
one can compute the color current. Inserting (|3.3|) into Eq. fl2.6|) , we get 

00 (— a )n+l 



3^) 



n=0 



nl 



+ (-1 



■ n+l 



d n f( Pai ) 
dp ai ■ ■ ■ dp an 



+ Tr[T a T Cl ---T Cn 



^/Tr[r a r Cl • 

d n fM 
dp ai ■■■dp 



dp ai ■ ■ ■ dp an 



(3.9) 
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The factor 2 in the above equation arises from the two helicities associated with every particle 
species. When the functions /, / and f g are determined, one can evaluate the momentum 
integral of Eq. (|3.9p, and then, after solving Eq. (|2.9[), one obtains the associated effective 
action. 



IV. PERTURB ATIVE SOLUTIONS 
A. General Considerations 



In Sec. Ill we have found exact solutions of the collisionless transport equations (|2.3| ) 
for translation invariant systems that obey the extra condition Eq. ( |3.7| ). In this section we 
treat the transport equations perturbatively and find solutions that are not constrained by 
the condition Eq. (|3.7|). 

We assume here that 

9 « 1 , (4-1) 

i.e. we deal with the weak coupling regime of the theory, and consider an expansion of the 
distribution function of the form 

Q = q(o) + q(D + q(2) + _ (4.2) 
where the 0-th term is a known function 

Q m (p,x) = f(p Qi ) . (4.3) 
The higher-order terms are determined by the equation 

p»D,QW(p,x) + |^ jiV(x), dQ(n d 2 P,X) } = • ( 44 ) 

Of course, the same treatments should be followed to study the antiquark and gluon distri- 
bution functions, but those are nearly identical. 

The iterative procedure based on Eq. ( |4.4|) does not correspond to a strict expansion 
in powers of g. Such an expansion would force us to split all covariant derivatives into a 
derivative and a commutator part, resulting in the breaking of the gauge covariance of every 
term in the perturbative series. We maintain the gauge covariance of every term in Eq. (|4.2| ) 
at the expense of reorganizing the perturbative expansion. 

It should be noticed that for / being the equilibrium distribution function, the first term 
in = 1) of the above perturbation series reproduces the hard thermal loops of QCD [fj, ^f]] . 
Let us also stress that while we are applying here a perturbative method to the transport 
equations in their matrix form, fully equivalent results can be obtained using the projected 
equations ( |2.5| ). Such an approach has been carried out by Bodeker and Laine to order g 2 for 
static quasi-equilibrium systems ||231 . We extend that analysis by pushing the perturbative 
procedure to higher orders in g. We do not attempt to solve the transport equations in full 
generality, as the solutions at every order then turn out to be highly non-local. We reduce 
our study to translation invariant systems when the solutions are much simplified. 
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B. From g to g 4 order 

In this subsection, we consider a translation invariant system that obeys Eq. ( |3.1|) for a 
fixed ccj. We also assume that the unperturbed distribution function depends only on p a . 
as in Eq. ( [4.3|) . The first-order correction to is obtained by solving the equation 

P"D,Q^ = -g p»F mi d -^- = - 9I fD,A ai . (4.5) 

Op ai OPa\ 

The solution reads (up to a function h, such that p^D^h = 0, that we will neglect through- 
out) 

QW(p, x ) = -gA ai (x)^^, (4.6) 

and it coincides with the first term of the ansatz (|3.3a| ) . 
The transport equation at second-order is 

P»D^ = £ {p»D,A ai ,A a2 } , (4.7) 



but it can be rewritten as 



M ^l^O A A a A m) ^. (4.8) 

Thus, the second-order solution is 

(p, X ) = g l Aai (x)A a2 (x) f /( ^ J , (4.9) 
2 op ai dp a2 

which again coincides with the respective term of the ansatz ( |3.3a| ). 
The third-order equation is 

P»D,Q® = -t- {p»D,A ai ,A a2 A a3 \ d3f iPa ^ , (4.10) 

i 1 J dp ai dp a2 dp r 



«3 



and its r.h.s. is not proportional to a covariant derivative of a third power of A ai . However, 
the anticommutator of the r.h.s. of the equation can be rewritten as follows: 

)(3) _ 9 3 /Oin (A a a \ i l|"r_.Mn A A 1 a \\ d 3 f(P° 



3! 



p»D„QW = -^[p»D lx (A ai A a2 A a3 ) + -\p>*D lx A ctl ,A a2 ],A a3 - J ™' . (4.11) 



dp ai dp a2 dp. 



Ct3 



The term with the commutator cannot be expressed as a total covariant derivative. Thus, 
the solution of the above equation contains two pieces: a term that is local in the gauge 
fields and a non-local part. Namely, 

Q (3 \x,p) = -£ A ai (x)A a2 (x)A a3 (x) d3f iPa ^ + Q®(x,p) (4.12) 

o. Op ai Op a2 Op a3 



with the non-local term obeying the equation 

,3 



(3) 



•L 

12 



[p^D^A ai , A a2 ], A, 



a.', 



d 3 f(Pa 



dp ai dp a2 dp, 



(4.13) 



One observes that when the condition Eq. fl3.7| ) is satisfied, = and we recover the so- 
lution of order g 3 of Eq. (|3.3a|) . But this term is non-zero under more general circumstances. 
We also notice that if 

\d^A a . } A aj ] = , (4.14) 

the non-local piece is proportional to g 4 . 

Let us solve Eq. (^4.13|) . Using the SU(N C ) algebra, we rewrite the commutator as 



,,/' n r)( 3 ) — — — f acd f deb T a (n^ F> A Y A c A b — 



= R a T° 



Consequently, = r a with 



p'd^qT = R a 



The solution can be expressed as 



>p- D< 



) ab R b (y) 



(4.15) 



(4.16) 



(4.17) 



where 1/p ■ D is the retarded Green function associated to the differential equation ( 4.16|) . 



The explicit form of 1/p ■ D can be found, for example, in Sec. II of [2S 
The fourth-order equation reads 
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9 4 /(Po 



OQ 



(3) 



12 



dp ai dp a2 dp a3 dp ai 24 



dp 



The first anticommutator in the r.h.s. of the equation can also be rewritten as a total 
covariant derivative plus an additional commutator. Thus, 



P^D^Q^ ^ — — ypf'D^ (A ai A a2 A a3 A a4 ) + [p^D^A^, A a2 ], A a3 A a4 

f 4! I ^ d Pu J ' 

d'fiPaJ 



d'f{p c 



dp ai dp a2 dp a3 dp 

(4.19) 



We see that the solution is of the form 



Q {A \x,p) = ^A ai (x)A a2 (x)A Q3 (x)A a4 (x) 



dp ai dp a2 dp a3 dp 

c 



(4.20) 



where represents the non-local contribution that vanishes if the condition Eq. (|3.7| ) is 
satisfied. We note that the local piece coincides with the term g A of the ansatz 
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In principle, we could solve the equation for and for the following terms of the 
perturbative expansion, finding that every term contains both local and non-local pieces in 
the gauge fields. However, we stop our analysis here as the structure of the non-local terms 
becomes more and more complex. 

Adding the contributions of quarks, antiquarks and gluons, we can obtain the color 
current. Up to the term n — 3, the color current still has a relatively simple form and is 
given by 



2E 



'-g) 



n+l 



n=0 



TV. 



A c a \(x) ■ ■ ■ AZ(x) JdPpf 



+ (-iy 



-i d"f(p ai 



dp ai ■ --dpo^. 
where the non-local term reads 



+ Tr[T a T Cl ...T Cn 



N f Ti[r a r Cl ---r ( 

d n f g ( Pai ) 



d n f( Po 



dp ai ■■■dp 



(4.21) 



dp ai ■ ■ ■ dp an 



+ 3aA x ) 



3a,n\( X ) 



9_ 
6 



+ 



j d pj d S{A^\y) 

d 3 f( Pai ) 



2N f Tr f r b 



dp ai dp a2 dp 



+ Tr (r' 



[P ' ^Alu -^O^]' ^ 



[p ■ DA ai , A a2 \, A a3 

d Z fg{Pa t 



d 3 f(p 



Cf.il 



dp ai dp a , 2 dp, 



dp ai dp a2 dp 

(4.22) 



In the following section we will get the effective action for static systems close to equilib- 
rium. In that case, it is easy to see that for n > 4 the local pieces do not contribute to the 
current or effective action, as either the momentum integrals (Appendix A) or the traces of 
generators in the adjoint representation (Appendix B) vanish. However, the same does not 
hold true for the non-local terms. 



V. EFFECTIVE ACTION OF QUASI-EQUILIBRIUM SYSTEMS 

In this section we find the effective action for systems close to thermal equilibrium. Then, 
the functions /, / and f g are of the Fermi-Dirac or Bose-Einstein form 



Jfb(E) 



Ifd(E) 



/be(£) 



1 



e /3(£-M) + 1 ' J*"\-~/ ©9(^+A*) + l' e /3E_l> 

where (3 = 1/T, and T is the temperature, and \i is the quark chemical potential. 



(5.1) 



A. Static Systems 

We first consider static systems satisfying the condition 

[D^A ,A ] = 0, fi = 0,1, 2, 3 



(5.2) 



The color current is given by Eq. ( |3.9| ) with = while the distribution functions are those 
in Eq. Q. 

After performing the momentum integral, one observes that j* = since 
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Air 







(5.3) 



where v % = p l /\p\ is the particle velocity, and the integral is performed over angular direc- 
tions of v. Thus, only j®, i.e. the color current density, is non- vanishing. In the case of static 
systems, one easily finds the Lagrangian density by integrating Eq. (|2.9|) . The fermionic 
contribution arises from 



C 



f 



7T 2 f- n 



-9) 



n+l 



n=o (" + l; 
while that of gluons is 



dEE 2 



d n f F p(E) 
dE n 



+ (-1) 



n+1 d n f FD (E) 



dE*- 



Tr[^ +1 (x)] , 



(5.4) 



Cn = Z 



7T 



oo 
n=0 



\n+l 



n + l 



dE E' 



d n f B E(E) 
dE n 



Tr[^ +1 (x)] 



(5.5) 



In the above expressions we have made use of the mass-shell condition, which gives E = 
Po = |p|, and we have put an infrared cut-off, T 3> A gT, which excludes the contribution 
of the soft modes to the integrals, as only the hard modes behave as quasiparticles. 
We note that the above Lagrangians can be written as 



C 



TN 



f 



Tr 



dEE 2 



ln(l + e-PiE-^-gA^ + ln(1 + e -f3(E+» +g A^ 

£ q = ±-Tr r dEE 2 ln(l - e -^-^ )) , 

7C 2 J A 



(5.6) 
(5.7) 



because 



d 



ln(l±e" 



dx y e x ± 1 

The Lagrangians ( |5.4j , |575| ) differ from those given in Eqs. ( |5.6| , |5.7|) by field independent 
terms corresponding to the free-quark and free-gluon contributions to the system's pressure. 

While it is not obvious from Eq. Q5.5D , Eq. (|5.7|) shows that the gluon action is not 
well defined without the infrared cut-off because for energies smaller than the eigenvalues 
of gAo, which are typically of order gT, the logarithm in Eq. (|5.7|) becomes a multi- valued 
complex function. Of course, there is an analogous problem with the Lagrangian Eq. ( |5.5| ). 
The difficulty one encounters is not surprising - the kinetic description breaks down for 
sufficiently soft modes. Then, the equilibrium distribution function of gluons of the form 
( p.3c|) becomes negative, loosing its probabilistic interpretation. 

We leave the evaluation of the energy integrals from Eqs. ( |5.4j ) and ( p.5| ) for the Appendix 
A, where we first compute the integrals with the cut-off. However, for all contributions the 
cut-off dependent terms appear to be regular and subdominant. Therefore, we take the 
limit A — > and thus drop all the cut-off dependent terms!. As a result, we find the field 
dependent fermionic contribution to the system's Lagrangian as 



§A matching procedure with the soft classical field theory allows one to eliminate the cut-off 
dependent terms but we will not pursue that procedure here. 
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% - '4(r- + + «-(^ + ^Al + ^ K + J^KAi . (5.8) 

We have kept here a linear term in the gauge potential that only survives in the Abelian limit 
where TrA = A . Taking into account that, as shown in Appendix B, the symmetric traces 
of an odd number of adjoint generators vanish, we find the following gluon contribution 

a 2 T 2 o A 

£ 9 = V^-^^" (5 - 9) 

It is also shown in Appendix B that the traces of the fields in the adjoint representation can 
be expressed through the fundamental representation traces. 

As already stressed, the integral ( |5.7| ) is ill-defined without the infrared cut-off. However, 
the integral with A = is regular after a Wick rotation to Euclidean space-time. For 
constant background fields A , and after the Wick rotation A — > we observe that 
Eqs. ( |5.6|) and ( |5.7|) with A = fully agree with the complete one-loop contribution to the 
effective potential for the phase of the Polyakov line 0,^,^]. It should be mentioned that 
a constant background field A can always be chosen in diagonal form , and thus satisfies 
the condition Eq. (|5.2|) . The Wick rotated integral ( f^7\ j is discussed in Appendix C. We 
note here that except for the terms analogous to those present in Eq. (|5.9| ) it provides non- 
analytic terms which are cubic in the eigenvalues of Aq . We also observe that Eqs. (|5.8| ) 
and ([5.9|) after the Wick rotation agree with the dimensionally reduced effective action . 

Now, let us briefly consider a static system that does not satisfy Eq. (|5.2|) . Using the 
solutions of Sec. [IV B| , we find that up to the order g 4 Eqs. ( p.8| ) and ( fT9| ) still hold, but at 
order g 4 and beyond there are corrections due to the non-local terms. The first non-local 
contribution appears at order g A from the non-local current Eq. ( f4.22p with a,i = 0. To get 
the corresponding term in the effective action, one should still solve Eq. ( |2.9| ). We have not 
found the solution, but we also see no reason why the term should vanish. Thus, we conclude 
that the presence of the non-local contribution signals a discrepancy with the dimensionally 



reduced effective action p7] at order g 4 , whenever the condition Eq. (p72|) is not satisfied 



B. Homogeneous Systems 

We consider here time-dependent homogeneous systems that obey the condition Eq. (|3.7| ) 
with oil = 1,2,3. For a system close to equilibrium, where the rotational symmetry is 
not broken, the solution of the transport equations can only depend on the modulus of 
the canonical three-momentum |p — gA(t)\. The color current in this situation is given 
by Eq. ([T9]) with f(p ai ),f(p ai ) and f g (p m ) replaced by /fd(|p|), 7fd(|p|), and /be(|p|), 
respectively. 

The color current is then given by Eq. ( |3.9| ), but now oti = 1, 2, 3. To get the expression 
for the color current, we first need to perform the momentum integral of Eq. (|3.9|) . The 
partial derivatives appearing in Eq. ( ^.9| ) can be explicitly evaluated. Using the chain rule, 
they can be expressed as partial derivatives of E = |p|. Therefore, 

df _ dE df d 2 f d 2 E df dE dE d 2 f 

dpi dpi dE ' dpidpj dpidpj dE dpi dpj dE 2 ' 
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d 3 f d 3 E df + fdE d 2 E ^ dE d 2 E ^ dE d 2 E \ d 2 f 



dpidpjdpk dpidpjdpkdE \dpidpjdpk dpj dpidpk dpk dpidpj J dE 2 
dE dE dE d 3 f 

where i, j, k = 1,2, 3. The partial derivatives of E can be written as functions of the particle 
velocities. Let us define 

d n E 
dp h ■ ■ ■ dpi 

An explicit evaluation of the first terms gives 



4r = - — hi • (5-io) 



v l v J 



L% = ~ (Lfav* + L$v> + Ifav') , etc. 



The angular integral of every term in the series of Eq. ( |3.9| ) can be easily calculated. 
Since the integrals of an even number of v 1 vanish we only need: 

*^,V = i**, (5.11) 



4?r 15 

The integrals over E are the same as those used in the previous subsection (see Appendix 
A). 

The evaluation of the momentum integrals gives, as expected, j® — 0. The space com- 
ponents of the color current are, of course, non-zero with the quark and gluon contributions 
given, respectively, as 

n 2 T 2 n 2 



+ JLjfTrfaA'A'A*) + Tr(T a A i A j A i ) + Tr(r a AMM J )) , (5.13) 



These currents arise from the following Lagrangian densities: 
and 

9 2 T* tA A , g 4 A A A A , 1 



fa = —^<T a A) + ^{Tr (T a A> A 1 A 1 ) + Tr(TJMM') + Tr (T a A i A i A^ ) j . (5.14) 
se currents arise from the 

jf = ~(y + ^(AjAA - -£- 2 (TriAjAjAAi) + ^{A^A)) , (5.15) 



C g = -^-Tr(^^) - \Tr(AjAjAiAi) + -Tr ( Ai A j A\ A j) ) . (5.16) 

We should point out that the terms proportional to g 2 in Eq. ( |5.15| ) and Eq. ( p.lfiQ 
agree with the hard thermal loop Lagrangians in the homogeneous limit |4[]. If the condition 
Eq. ( |3.7| ) is not satisfied, the color current, and thus the associated effective Lagrangians, 
are corrected at order g 4 and beyond by the addition of the non-local terms, exactly as 
happened for the static systems. However, we will not write down these explicit terms here. 
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VI. CONCLUSIONS 



Our results show the efficiency of transport theory in describing the quark-gluon plasma 
at soft scales. We have shown how the solutions of the collisionless transport equations for 
the static systems close to equilibrium reproduce the one-loop effective potential for the 
phase of the Polyakov line. Up to now, the results of the transport theory and quantum 
field theoretical computations have been known to agree only for the lower-dimensional 
operators, but our computations indicate that the agreement extends to the full one-loop 
effective action. We find a complete match of the transport results with those of the one-loop 
effective potential in the presence of a constant background field. For non-constant static 
background fields, transport theory predicts the appearance of non-local operators in the 
effective action, starting at order g A and beyond. This result then suggests a discrepancy 
with the dimensionally reduced effective theories ||27|| , where these non-local operators are 
not present. 

We have limited our analysis to translation invariant systems, when the solutions of the 
transport equations are local in the lower orders of the perturbative expansion. It would be 
desirable to solve the equations in full generality. However, the solutions are then complex 
non-local functions of the gauge fields, and their structure beyond order g 2 is not particularly 
enlightening. It is presumably more promising to explore the combined set of equations ( |2.3| ) 
and (p.8| ) numerically, thus allowing for a non-perturbative study of dynamical phenomena 
at soft scales, even beyond the hard thermal loop approximation. 
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APPENDIX A: MOMENTUM INTEGRALS 
1. Bosons 

The bosonic integrals to be evaluated are of the form: 

I b (A) = T 3 -" Hdxx 2 d -M^L , (Al) 
n 7a dx n 

where A = A/T. Expanding the Bose-Einstein distribution as 

I oo 

/be(x) = = E e~ mX 

e - 1 - m=l 
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and interchanging the order of summation and integration, the integral that has to performed 
reduces to 

' 2 2A A 2 \ 



dxx e 



-mA 



ni 



Therefore, 



I»(A) = 2(-l) n T s ~ n Li 3 _„(e- A ) + ALi 2 _ n (e" 



~~2 



where the solution is expressed in terms of the Euler polylogarithm function 



oo m 

U s {z) ^ V — . 

m=l 



The function Li^e A ) can be expanded in powers of A as [ |5D|j5T| ] 

oo 



LW(e" A ) = E(-i: 

k=0 



for N < 1 and 
U N (e~ A ) 



E 

k=0, k^N-1 

for N > 1; the zeta function is defined as 



(lnA-^.x) 



(A2) 



(A3) 



(A4) 



(A5) 



(A6) 



c(«) = E 



m=l 



Res > 1 . (A7) 
|, A6) with the zeta function 



and H N = 1 + \ + \ + • • • + ^. 

While the series (|A4]) is convergent for \z\ < 1, the series ( 
defined by Eq. (|A7|) seem to be divergent because the zeta argument is repeatedly equal or 
smaller than 1. This happens due to the 'illegal' interchanging of the two summations. As is 
well known, the problem is resolved by means of the zeta function regularization procedure 
P^,|3l| , where the analytic continuation of ((s) instead of the definition ( |A7D is used. Then, 
C(l) remains truly divergent while 



C(Q) = -\ , 



C(l-2*) 



B 



2k 



2k 



k = 1,2,3, 
k — X . 2, 3 . 



a-2k) = o , 

where Bi are the Bernoulli polynomials. 

Using Eqs. ( |A5| , |A6| ) one finds that when A — > 

U N (e~ A ) = ((N) + 0(A) (A8) 

for N 1. Since Lii(z) = — ln(l — z), Li^e - ^ diverges as — InA. Therefore, when A — > 

/^2(-l)"T 3 -"C(3-n) (A9) 

for n ^ 2. The n— even contributions to the effective action vanish anyway because the 
respective trace of the A-fields equals zero, see Appendix B. 
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2. Fermions 



The fermionic integrals of interest are 

Il{a, A) = T 3 ~ n (4(a, A) + 4(-a, A 

where a = j3[i and 

d n f fd(% ~ a) 



Jl(< 



dxx z 



dx n 



(A10) 



(All) 



Expanding the Fermi-Dirac distribution as 

1 



/ FD (x - a) = = 

e x + 1 ~! 



m-lg-m(i-a) 



and interchanging the order of summation and integration, we obtain after performing the 
integral 



(A12) 



J^a, A) = 2(-l)"- 1 ^Li 3 - n (-e a - A ) + ALi 2 _ n (-e a - A ) + ^Li 1 _„(-e a - A ) 
Using the formula p2) 



where rj(s) is the alternating zeta function defined as 

oo ( j\m— 1 



m=l 



(l-2 1 -)C(a) 



Res > 1 



we find that 



00 a 1 



Jf{a, A) = 2(-l)« £ - 7/(3 - n - I) + 0(A) , 

when A — > 0. 

Inserting the series ( |A15| ) into Eq. (|A10 ) one gets 

°° n l r 

//(a, 0) = 2T 3 - 2 7T ^(3 - n - [(-l) n - (-1)' 



(A13) 



(A14) 



(A15) 



(A16) 



As seen in Eq. (|A16 ), the argument of rj is always an even number for non- vanishing 
terms. Since rj(—2k) = for k = 1,2,... the series in Eq. QA16 ) terminates. Then, one 
observes that l£(a, 0) = for n > 4, while 
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Il{a, 0) = T 3 ( 4a rj(2) + |a\( 0)) = T 3 ^ a + I a 3 " 



//(a, 0) = T 2 (-4 77(2) - 2a 2 77(0)) = —T 

j/(a,0) = 4aTr7(0) = 2aT , 
/ 3 / (a,0) = -4r/(0) = -2, 



3 



(A17) 



where 77(2) = vr 2 /12 and 77(0) = 1/2. 



APPENDIX B: ADJOINT REPRESENTATION TRACES 



To compute the effective action of Eq. ( |5.8|) and Eq. ( |5.9|) one needs to evaluate traces in 
the fundamental and adjoint representations. We present here some useful formulas which 
relate the traces of the fundamental and adjoint representations. 

First, it is easy to prove that the total symmetric traces of an odd number of adjoint 
generators vanish. In order to prove this, note that 

Tr [A?] = Tr[(^f] , (Bl) 

where the superscript T denotes transposition. The adjoint representation of SU(N C ) is 
real, and the generators obey Tj = — T a . Therefore, 

Tr[A£] = (-l) m Tr[^™] . (B2) 

Thus, the trace vanishes for odd m. 
Two other useful formulas are: 

Tr„4 2 = 2N c TtA 2 , (B3) 
Tr^ = 6(TrA 2 ) 2 + 2N c TrA* . (B4) 

For N c = 2 and N c = 3, one also has 

TrAt = i(TrA 2 ) 2 . (B5) 



APPENDIX C: THE ONE-LOOP BOSONIC ACTION 



In this appendix we show how the solutions of the transport equations allow one to 
reproduce the results of P, |25| , |26| . We consider only the bosonic case for the SU{2) gauge 
group. The fermionic one can be treated similarly. Due to a global color rotation a constant 
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background field Aq can be always chosen in the diagonal form Aq = S a3 C/g where C is a 
real constant. The respective gauge field in the adjoint representation then reads 



A 



Now, we plug this expression in Eq. 
Taking the trace ( Tr^ m = 2(C/g) 2m , 
resumming reads 





and expand the logarithm and the exponential. 



TrA 



2m+l 



0), we end up with a series which after 



C =1 

9 7T 2 



dE E 2 



Rotating to Euclidean time Ac 



T 



7T Z 



dE E 2 



ln(l - e 



-/3(E-C)^ 



+ ln(l - e 



-0{E+cy 



iA E , and thus C 



iC E , we find 



ln(l 



-/3(£-iC fe )> 



+ ln(l 



-P(E+iC E )\ 



(CI) 



(C2) 



25] . The integral is evaluated in 



which matches the result found in Eq. (23) of 
and non-analytic cubic terms in C E are found. When rotated back to the real time these 
terms provide imaginary contributions to the effective action. Its presence in the Minkowski 
effective action could have been anticipated from Eq. (|CTD, as for soft energies E < C the 
logarithm becomes a multi-valued complex function. The real time integral is not well- 
defined and an additional prescription, such as going to Euclidean time, is necessary in 
order to evaluate it. With the infrared cut-off A > C, these problems are absent, and the 
Minkowski effective action is real. 
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